EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2012-164 
2012/07/04 



CMS-EXO-1 1-067 



Search for new physics with long-lived particles decaying 
to photons and missing energy in pp collisions at 

= 7TeV 



The CMS CoUaboratioif] 



Abstract 



A search is performed for long-lived neutral particles decaying into a photon and in- 
visible particles. An example of such a signature is the decay of the lightest neutralino 
with nonzero lifetime into a gravitino and a photon in gauge-mediated supersymme- 
try with the neutralino as the next-to-lightest supersymmetric particle and the grav- 
itino as the lightest. The search uses events containing photons, missing transverse 
energy and jets. The impact parameter of the photon relative to the beam-beam col- 
lision point can be reconstructed using converted photons. The method is sensitive 
to lifetimes of the order of 0.1 to 1 ns. The data sample corresponds to an integrated 
luminosity of 2.23 fb^^ in pp collisions at i/s = 7TeV, recorded in the first part of 2011 
by the CMS experiment at the LHC. Cross-section limits are presented on pair pro- 
duction for such particles, each of which decays into a photon and invisible particles. 
The observed 95% confidence level limits vary between 0.11 and 0.21 pb, depending 
on the neutral particle lifetime. 
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1 Introduction 

New heavy particles with long lifetimes are predicted in many models of physics beyond 
the standard model (SM), including models with gauge-mediated supersymmetry breaking 
(GMSB) |T-'9l and certain classes of hidden valley models [10|. Such particles may be neutral 
and decay into photons and invisible particles, which escape detection. Their lifetime is es- 
sentially a free parameter of the models, which motivates searches sensitive to a wide range 
of potential lifetimes. For sufficiently long lifetimes, of order of 1 to 10 ns, measurement of 
the shower direction by the DO experiment [TT] or time-of-flight by the CDF experiment ||T2| 
with an electromagnetic calorimeter can be used to identify such decays, and limits have been 
obtained by using such techniques at the Tevatron. 

In this paper, we present a new analysis method that is sensitive to neutral particles with life- 
times of the order of 0.1 to Ins. In a fraction of events, photons undergo conversion into 
e+e^ pairs in the material of the LHC beam pipe or the CMS charged-particle tracking system, 
providing a clear experimental signature. The tracks of the electrons can be precisely recon- 
structed and used to calculate the photon trajectory. In particular, the photon impact parameter 
with respect to the primary vertex can be determined. 

We search for diphoton events with at least one converted photon having a significant im- 
pact parameter produced in association with missing transverse energy. Models with gauge- 
mediated supersymmetry breaking, e.g., general gauge-mediation (GGM), are used to demon- 
strate the kinematics of the production and decay. Figure [l] shows one of the production and 
decay diagrams. Assuming that R-parity is conserved [13], SUSY particles are produced in 
pairs and decay into SM particles and the lightest neutralino (^). The neutralino decays into 
a photon and a gravitino (G), the lightest SUSY particle in this model, which escapes the de- 
tector, leading to apparent missing transverse energy £^'**. Moreover, in many models, the ^ 
is produced in association with high transverse momentum (px) jets. We consider neutralino 
proper decay lengths between ct = 2 and 25 cm, corresponding to a lifetime t on the order of 
0.1 to Ins. 




Figure 1: Example diagram of Xi pair production and ~^ 7 + G decay. 

The data sample was collected in 2011 by the Compact Muon Solenoid (CMS) experiment at 
the LHC from proton-proton collisions at a center-of-mass energy of 7 TeV, corresponding to 
an integrated luminosity of 2.23 ± 0.05 fb^^. A detailed description of the CMS detector can 
be found elsewhere [14J. The detector's central feature is a superconducting solenoid pro- 
viding a 3.8 T axial magnetic field along the beam direction. Charged particle trajectories are 
measured by a silicon pixel and strip tracker system, covering < ^ < 27r in azimuth and 
\rj\ < 2.5, where pseudorapidity is defined as i] = — In tan 0/2, and 6 is the polar angle with 
respect to the counterclockwise beam direction. The amount of material crossed by a particle 
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2 Event sample and selections 



traversing the tracker volume from the interaction region varies from 0.4 to more than 2 ra- 
diation lengths, depending on the particle direction. A lead-tungstate crystal electromagnetic 
calorimeter (ECAL) and a brass-scintillator hadron calorimeter (HCAL) surround the tracker 
volume. For the barrel and endcap calorimeters ( | ?/ 1 < 3.0), the modules are arranged in projec- 
tive towers. Muons are measured with gaseous detectors embedded in the steel return yoke of 
the magnet. The detector is nearly hermetic with calorimeter coverage up to | ?/ 1 = 5.0, allowing 
reliable measurement of E^^^. 

The analysis strategy consists of selecting events with a diphoton final state with at least one 
converted photon, and then examining the impact parameter of each single photon for a dis- 
placed photon signal. Missing transverse energy and at least two jets are also required. The 
background is estimated using a low E^^^^ control sample. Upper limits on the cross section for 
pair production of neutral particles, each of which decays into a photon and invisible particles, 
are then computed as a function of the neutral particle's lifetime. 

2 Event sample and selections 

The data were recorded using the CMS two-level trigger system. This analysis selects events 
with at least two photons. A diphoton trigger is required with ECAL transverse energy thresh- 
olds Ej set to values increasing from 32 (22) to 40 GeV (28 GeV) for the leading (sub-leading) 
photon, as the instantaneous luminosity rose over the course of the data-taking period. To en- 
sure being on the plateau of the trigger efficiency, the offline analysis selects events with at least 
two photons with Ej > 45 GeV (30 GeV) for the leading (sub-leading) photon in the event. The 
data sample is used both for the selection of signal candidates and for control samples used for 
background estimation. 

The PYTHIA 6.4 event generator |[T5l | is used to simulate SUSY signal events. In particular, we 
generate SUSY GMSB signal datasets in the SPSS benchmark model fl^ with A around 100 TeV. 
Every event in this sample has two neutralinos, each having a mass of 140 GeV. The neutralinos 
in each dataset have a specific mean lifetime. The signal selection efficiency can depend on the 
^ mass, but for masses around 100 GeV this dependence has a negligible effect in this analysis. 
Each of the neutralinos decays to a photon and a light gravitino. The CMS detector response is 
fully simulated using Geant4 [17J. 

Because of the high luminosity during the run, many events had multiple pp interactions lead- 
ing to a distribution of primary vertices, called "pile-up". The pile-up conditions at high lu- 
minosity affect the event selections, in particular conversion reconstruction efficiencies. The 
generated pile-up distribution in the Monte Carlo samples was reweighted to reproduce the 
2011 data-taking conditions. 

The unconverted and converted photon candidates are reconstructed from clusters of energy 
in the ECAL. At least one photon candidate with transverse energy Ej > 45 GeV and re- 
constructed in the ECAL barrel region is required. The transverse distribution of energy in 
the associated cluster of ECAL crystals must be consistent with that expected from a pho- 
ton, and the energy detected in the HCAL behind the photon shower cannot exceed 5% of 
the ECAL energy. To suppress hadronic jets misreconstructed as photon candidates, we re- 
quire the latter to be isolated from other particles in the tracker, ECAL, and HCAL. A cone 
of AR = (Arj)^ + {A(p)^ = 0.4 is constructed around the candidate's direction, and the 
scalar sums of transverse energies of tracks and calorimeter deposits within this cone are 
determined, after excluding the contribution from the candidate itself. The isolation sums 
are required to be E Cecal < 0.006 ■ Ej + 4.2 GeV, EErcal < 0.0025 ■ Ej + 2.2 GeV, and 
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^ptocks < 0.001 • Et + 2.0 GeV, for the ECAL, HCAL, and tracker, respectively, where Ej is 
the transverse energy of the photon candidate. To reject electrons, we require that photon can- 
didates do not have associated hit patterns in the pixel detector consistent with an electron 
track and not be part of a reconstructed conversion track pair. 

Photon-like signals that satisfy the same requirements as photons except for isolation are de- 
fined as "fake photons". These objects, predominantly jets containing energetic neutral parti- 
cles such as tt'^s, are used for the background estimation. 

Jets are reconstructed from energy deposits in the calorimeters using the anti-A;x clustering al- 
gorithm |T8| . The energies of these jets are corrected using the pj of the charged tracks re- 
constructed in the tracker |19|. At least two jets are required with \tj\ < 2.6 and leading jet 
transverse momentum pxi > 80 GeV and sub-leading jet transverse momentum px2 > 50 GeV. 

Missing transverse energy is calculated from calorimeter energy deposits. It is corrected us- 
ing tracking information for energy missed owing to incomplete calorimeter measurements of 
muons and charged hadron energies, especially for soft tracks that do not reach the calorime- 
ters Ii20i| . The E™^^ threshold of 30 GeV discriminates between signal and background. 



3 Photon conversion and the photon impact parameter 

The CMS silicon tracker provides robust and precise reconstruction of charged-particle mo- 
menta in the high-occupancy environment of LHC collisions. Photon conversions in the tracker 
material can be used to obtain the photon direction. By extrapolation along the momentum di- 
rection from the conversion vertex back to the beam line, we can calculate the impact parameter 
of the displaced photons, as shown in Figure |2] If the has a nonzero lifetime, the decay pho- 
ton can originate from a displaced secondary vertex rather than close to the primary vertex. 
The momentum direction determined from conversion tracks for this photon will point away 
from the primary vertex, giving a nonzero impact parameter (IP), which could be a signature 
for a long-lived signal. The conversion reconstruction efficiency is approximately 5%. 



G 




Primary vertex 




Primary vertex 



Figure 2: Left: Illustration of ^ — )■ 7 + G decay. Right: In x-y plane, a neutralino travels 
with a finite lifetime and decays into 7 + G. The diagram shows the photon converting into an 
e+e^ pair, and the subsequent reconstruction of the impact parameter. 



The transverse impact parameter dxy is the distance of closest approach of the photon trajectory 
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3 Photon conversion and the photon impact parameter 



to the beam line in the transverse plane to the beam. The longitudinal impact parameter dz is 
the distance of closest approach of the photon trajectory to the primary vertex projected on to 
the z direction. The photon trajectory is defined as a straight line from the conversion vertex 
along the conversion momentum. The impact parameters dxY and are expressed in terms 
of the vector L = (Lx, Ly, Lz) between the conversion point and the primary vertex, and the 
angle a of the conversion momentum vector p projected on to the x-y plane as follows: 



d-xY = —Lx ■ sin a + Ly ■ cos oc 

Lx ■ px + ■ Py Pz (1) 



d 



z 



Pi Pi 



The momentum p is calculated as the vector sum of the e+e pair momenta at the conversion 
vertex. 

The e^e^ tracks at the conversion vertex are produced with a small opening angle. As a kine- 
matic constraint on their Q and (p angular directions, the momenta of the tracks are required to 
be parallel. This constraint is applied to select conversion tracks. Similarly, the conversion ver- 
tex is fitted using kinematic constraints. Three algorithms are used to reconstruct conversions: 
the "tracker-only" [21, 22J, the "ECAL-seeded" and the "Gaussian sum filter" (GSF) EH. 
The tracker-only algorithm selects track pairs from all reconstructed tracks, with the kinematic 
constraint, and fits the conversion vertex; the ECAL-seeded algorithm takes the energy deposits 
in the ECAL as seeds and then extrapolates back to the tracker to fit the conversion tracks and 
vertices; the GSF electron algorithm follows a procedure similar to the tracker-only algorithm 
but uses electrons identified by the GSF algorithm to find the track pairs. The conversions from 
these three algorithms are merged and the duplicates are removed. Among the reconstructed 
conversion tracks, the one with the largest number of hits is kept if it passes the necessary 
quality cut. To be selected, the tracks of two opposite-signed conversion particles both need at 
least five valid hits; in addition, the conversion vertex requires a valid fit with probability 
> 5 X 10-4. 

In high luminosity conditions, multiple collisions occur within a single bunch crossing, pro- 
ducing a number of primary vertices. The true primary vertex can have a large deviation from 
the interaction point in the longitudinal direction but a much smaller uncertainty in the trans- 
verse direction. To be robust against pile-up conditions, the transverse impact parameter [dxr) 
with respect to the transverse position of the beam line is used in this analysis. To illustrate 
the dxy signature, the distribution of its absolute value in data for events with E™^^ > 30 GeV 
is compared with a Monte Carlo sample with ^ proper decay length ct = 5 cm, as shown in 
Figure 3j The requirement that I^xyI be larger than 0.6 cm then defines the signal region. This 
cut on I dxy I is derived by optimizing the expected limits on the production cross section. As 
seen in Figure |3| one event passes all selection criteria. The Monte Carlo sample is normalized 
to the integrated luminosity, and on average 8.3 events from the SPSS benchmark model are 
expected to pass all selections. 

The effect of the event selection is illustrated in Table [l] for the ct = 5 cm sample. Of the 45 057 
simulated events, 711 events survive the cuts, for an overall event selection efficiency of 1.58%. 
Efficiencies for four neutralino lifetimes are given in Table |2] On average, the signal efficiency 
is about 1% to 2% after the trigger requirements and selection criteria are applied. Beyond ct 
= 25 cm, the signal efficiency decreases due to the limitations of the tracking algorithms for 
reconstructing highly displaced tracks. 
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Figure 3: The |dxy| distribution for data with Ef^^ > 30 GeV compared with signal simulation 
for CT = 5 cm normalized to the integrated luminosity of the data. The arrow indicates the | dxr \ 
requirement of the signal region. 



Table 1: Signal selection cut flow for ct = 5 cm. 



Selection 


Events in Monte Carlo 


Total 


45057 


Diphoton trigger 


39988 


Leading (sub-leading) photon 
£t > 45 (30) GeV 


37398 


Any ECAL barrel photon Ej > 45 GeV 
and photon identification 


27766 


Jets Pti > 80 GeV and pj2 > 50 GeV 


26229 


Conversion selection 


1602 


E^'^^ > 30 GeV 


1542 


cfxY > 0.6 cm 


711 



Table 2: Event selection efficiency vs. proper decay length. 



CT [cm] 


2 


5 


10 


25 


Efficiency 


0.92% 


1.58% 


1.80% 


1.39% 


Statistical errors 


0.05% 


0.06% 


0.06% 


0.06% 



4 Background determination 

The background contributions to the sample of events with photon-plus-jets with missing 
transverse energy arise mainly from QCD events with real photons, and QCD multijet events, 
due to the effect of finite detector resolution on the measurement of £^'®®. In QCD events with 
real photons, there are typically one real photon and one jet that is misidentified as a pho- 
ton in the final selection. The rate of QCD multijet events is sufficiently large to have a finite 
probability for jets to be misidentified as photons. There are additional sources of background 
with intrinsic missing transverse energy, such as the W+jets process, but these are heavily sup- 
pressed by the Ej cut for the photon and make a negligible contribution to the final selection. 
The ^ — ?■ 7 + G decay has two signatures: £^^^ from the unseen G and a large impact pa- 
rameter from the displaced photons. The strategy for determining the background is to use 
control samples that are kinematically similar to the candidate sample while having no real 
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5 Systematic uncertainties 



E^®*. The data with E^^^ < 20 GeV are selected as such a control sample. The background 
is then estimated as the number of events with £™^®® < 20 GeV and satisfying all the other 
selection criteria. 

To ensure that the E^'*® < 20 GeV region of data gives a correct estimation of the background, 
the \dxY \ distributions of fake photons, defined previously, and isolated photons are compared. 
First, the \dxY \ distribution of fake photons for events with E^'®^ < 20 GeV (background region) 
is normalized by the total number of conversions and compared with the [dxrl distribution of 
fake photons in the region where E^'*® > 30 GeV (signal region), as shown in Figure |4| top). 
The I^xyI distribution in the E^^^ < 20 GeV (background) region is reweighted by the con- 
version vertex probability of the E^®* < 20 GeV (signal) region. This takes into account 
that the distributions of the probability are not the same in the background and signal re- 
gion. Second, the |dxy| distributions of isolated photons and fake photons in the low E^^^^ 
region are compared by normalizing to the total number of conversions, also reweighted by 
the conversion vertex x^ probability, as shown in Figure [ij^bottom). The [dxrl distributions of 
fake photons and isolated photons agree, and thus we conclude that the | dxY \ distribution for 
events with E^^^^ < 20 GeV gives a good description of the background. The residual differ- 
ences in the fake and isolated photon distributions in the low E™*® control sample are consistent 
with the expected ervhancement of jets misreconstructed as photons in the fake photon sample. 
We therefore take the isolated photon sample in the low E^^*^ control sample as the basis for 
the background estimate in the signal region, after we apply conversion vertex x^ probability 
reweighting. 

The isolated photon {dxrl distributions for the signal and background regions are shown in 
Figure |5] They are normalized to the same number of conversions and reweighted by the 
conversion vertex probability to predict the number of background events. A requirement 
that Icfxrl > 0.6 cm is applied to the background region E^^^ < 20 GeV, which gives a total 
background of 0.78^'^q events. As a cross check, Monte Carlo background samples are fitted 
with a general functional form to extrapolate the expected tail of events for isolated photons 
beyond the \dxY\ > 0.6cm cut. Applying this general functional form to the background sam- 
ple derived from data and integrating the expected rate above the |dxy| > 0.6 cm cut gives a 
background estimate consistent with the number quoted above. 



5 Systematic uncertainties 

Table |3] summarizes the systematic uncertainties affecting this analysis. To determine the un- 
certainty in the conversion reconstruction efficiency, we compared the predicted and observed 
numbers of conversions for Z — > f//^7 events. For photons with Ej > 20 GeV, a scale factor of 
0.87 ± 0.08 is derived. High Ej photons have a reconstruction efficiency of (7.0 ± 0.5)% in the 
Drell-Yan to dimuon Monte Carlo samples, and the corresponding high Ej photon efficiency 
in the GMSB Monte Carlo samples is (6.0 ± 0.7)%, yielding a Z ^ mij-GMSB scale factor of 
1.16 ± 0.08. Considering the correlations between the conversion reconstruction efficiency and 
the high- Ex dependence, the uncertainty of conversion reconstruction is determined to be 21% 
for the signal region. 

The dxY resolution is studied by comparing the dxY resolution for converted photons from 
Z — )■ jifi'Y events between Monte Carlo and data. The data and Monte Carlo samples are in 
agreement and a dxY resolution of 0.06 cm extracted from data gives less than 0.5% uncertainty 
for the signal selection. The resolution uncertainty includes tails in the fitted resolution func- 
tion. 
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Figure 4: Top: Fake photon \dxY\ distributions in the background (£^'®^ < 20 GeV) and signal 
(£^iss y 30 GeV) regions. Bottom: Isolated and fake photon \dxY\ distributions in the back- 
ground region (£^'®® < 20 GeV). 
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Figure 5: \dxY\ distributions for isolated photons: background region (£^'^^ < 20 GeV) com- 
pared to the signal region (£^^* > 30 GeV). 
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7 Summary 



Other sources of uncertainty include the integrated luminosity (2.2%), the jet pj/E^^^ selection 
requirement, and photon identification. The jet pj/E^^^ resolution contributes less than 0.5% 
to the total uncertainty llT9l 1201 . The photon identification efficiency is studied in [25] and 
includes the uncertainties from pile-up (2.5%), the photon data/MC scale measured in a Z — > ee 
data sample (2.6%), and the photon-electron identification difference as studied in MC samples 
(0.5%). Other sources of uncertainty give negligible contributions to the systematic uncertainty. 
The total systematic uncertainty is determined to be 25% by considering all the correlations 
among pile-up, photon identification data/MC scale, and conversion reconstruction efficiency. 



Table 3: Summary of systematic uncertainties. 



Systematics 


Uncertainty (%) 


Integrated luminosity 


2.2 


Jet Pt/£^'*® energy scale 


< 0.5 


Pile-up 


2.5 


Photon identification data/MC scale 


2.6 


Photon-electron difference 


0.5 


Conversion reconstruction efficiency 


21 


Photon dxY resolution 


< 0.5 


Total 


25 



6 Results and interpretation 

One event with I^xyI = 0.74 cm and £^'^^ = 44.9 GeV satisfying all the other selection criteria 
is observed. The estimated background is O.ZS^Q ^g events. We determine the upper limits 
for the cross section for pair production of neutral particles, each of which decays into one 
photon and invisible particles. A CLs limit setting method p6l |27| is employed using log- 
normal uncertainties for the total background rate to incorporate the uncertainties in the total 
background rate, integrated luminosity, and total acceptance times efficiency. The observed 
95% confidence level limits vary between 0.11 and 0.21 pb, depending on the neutral particle 
proper decay length (Table |4] and Figure |6|. 

Table 4: 95% confidence level (CL) upper limits on the cross section for pair production of 
neutral particles, each of which decays into a photon and invisible particles, as a function of 
the neutral particle proper decay length. 



CT (cm) 


2 5 10 25 


a (pb) 95% CL 


0.21 0.12 0.11 0.14 



7 Summary 

We have introduced a novel method using the photon conversion impact parameter to search 
for new physics involving long-lived particles decaying into photons. The high resolution of 
the impact parameter reconstruction enabled our analysis to probe for new physics to much 
lower values of £^^*^ compared to other LHC searches [25|. The search was performed using 
the final state of photons, jets, and missing transverse energy. Cross-section limits on pair 
production for such particles, each of which decays into a photon and invisible particles, were 
set as a function of the long-lived particle's lifetime. These upper limits are applicable to a 
general class of new physics processes with diphotons in the final state. 
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Figure 6: 95% confidence level upper limits on the pair production cross section for neutral par- 
ticles, each of which decays into a photon and invisible particles, as a function of the neutral 
particle proper decay length. The observed values as a function of mass are shown by the solid 
line. The dashed line indicates the expected median of results for the background-only hypoth- 
esis, while the green (dark) and yellow (light) bands indicate the ranges that are expected to 
contain 68% and 95% of all observed excursions from the median, respectively. 
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